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1 Introduction 

The measured value of = r^ij/Fhad in decays is drawing considerable at- 
tention in the physics community. On the one hand, according to all four LEP 
experiments, it deviates significantly from the Standard Model prediction, a result 
that has remained quite stable in the last few years On the other hand, it is a 

considerably difficult measurement, and it is certainly disturbing that the only devi- 
ation from the Standard Model is found in a quantity that is affected by complicated 
strong-interaction effects. In the present paper we examine one possible source of 
systematic errors that may affect the measurement of R^, whose origin is purely dy- 
namical. In several experimental techniques for the measurement of R^, the tagging 
efficiency is extracted from data by comparing the sample of events in which only one 
b has been tagged, with the one in which both 6's are observed. If the production 
characteristics of the b and the b were completely uncorrelated, this method would 
yield the exact answer with no need for corrections. Of course, other correlations of an 
experimental nature should be properly accounted for, but their discussion is outside 
the scope of the present theoretical paper. Standard QCD gluon emission generates 
a correlation of the quark-antiquark momenta of the order of as- Other dynamical 
effects, such as the production of heavy quark pairs via a gluon-splitting mechanism, 
may affect the measurement. However, they are to some extent better understood. 
In the first place they tend to give soft heavy quarks, and they are therefore easily 
eliminated. Furthermore, experimental studies of these production mechanisms have 
begun to appear 

The aim of the present work is to study the dynamical correlation of the heavy 
flavoured quark-antiquark pair, with particular emphasis upon its impact on the de- 
termination of Rb- We will use leading-order QCD formulae throughout, emphasizing 
the cases in which a higher-order calculation would be desirable and possible. We 
found that a particularly relevant quantity is the average momentum correlation r 
defined as 

_ {X1X2) - {xf 
{xf ' 

where Xi^2) are the Feynman x of the produced B{B) mesons, Xi(2) = Pb{b)/p^b^^^\ 
and (x) = (a;i(2)). In fact, we will show that under the assumption that the detection 
efficiency for a B meson is nearly proportional to its momentum, the quantity r 
represents the relative correction to R^ due to momentum correlations. The quantity 
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r has a well-behaved perturbative expansion, starting at order as, in spite of the 
fact that the perturbative expansions of {X1X2) and (x) are not well-behaved, since 
their coefficients are enhanced by logarithms of the total annihilation energy over the 
heavy-quark mass. Our perturbative calculation for r gives a value of the order of 
1%. Power-suppressed corrections may still affect this quantity to some extent. If 
power corrections of the order of 1/Q are present, they may turn out to be of the 
same order of the perturbative value at LEP energy. Monte Carlo models ||^ seem to 
support this possibility p!0| . 

In general, the QCD formalism allows us to compute the double differential cross 
section for the inclusive production of the heavy-flavour pair, using as input the 
value of as and the fragmentation function of heavy-flavoured mesons. It is easy to 
compute such quantities in leading order. The leading order result may be improved 
on one hand by including next to leading corrections, and by resumming terms that 
are enhanced logarithmically near the end of the phase space (i.e. for xi(2)— At 
present these calculations are not available. 

Shower Monte Carlo models include to some extent perturbative effects of leading 
order, and effects that are enhanced near the end of the phase space. At the parton 
level, they should therefore provide a description that is compatible with the leading- 
order calculations reported in the present paper. Some recent studies 0] support this 
conclusion. On the other hand, fixed order matrix elements Monte Carlo programs 
should be used cautiously when computing correlations, because of problems arising 
in the truncation of the perturbative expansion. 

The rest of the paper is organized as follows. In Section Q we discuss the average 
momentum correlation, and how it affects the determination of Rf,. In Section |^ we 
describe the leading-order calculation of the double inclusive cross section for heavy 
flavour production. Finally, we give our conclusions in Section ^ 

2 The average momentum correlation 

We begin by considering the simple case in which the efficiency for tagging a B 
meson is a linear function of its momentum. This simplifying assumption allows us to 
make very precise statements about the correlation. Furthermore, it is not extremely 
far from reality, in the sense that the experimental tagging efficiency is often a growing 



Figure 1: Some of the diagrams contributing to the process Z^bb. Virtual 
diagrams are not shown. 



function of the B momentum. Let us therefore assume 

€b = Cxi, ei = Cx2, (2.1) 

where eb is the tagging efficiency. We assume that in our ideal detector the detection 
efficiency is not influenced by the presence of another tag. Therefore the total number 
of tags is then given by 

^ = R'bC i{x^) + {X2)) = R^C {x) , (2.2) 

where X2 is the Feynman x of the produced B meson, and with (x) we indicate the 
common value of their average. Observe that the in the above equation A^^i is the 
number of tags, given by the number of events with one tag plus twice the number 
of events with two tags. On the other hand, the number of events with two tags will 
be given by 

^ = R'bC = R'b {C {x)f X (1 + r), (2.3) 

where r is deflned in eq. (|1.1|) . We can then extract R^ 

Rl = — ^ X (1 + r) . (2.4) 

^ 4:NN2 ^ ' ^ ' 

The quantity r cannot be measured; one therefore has to compute it in order to 
determine 

The basic QCD processes for the production of a bb pair are depicted in flg. |T[ Let 
us focus on the process up to the order a^. The flnal state is fully deflned in this 
case by the variables xi = 2eb/E and X2 = 2eb/E, where E is the total centre-of-mass 
energy. The gluon momentum fraction X3 is given hj 2 — xi — X2- We can write the 
differential cross section, normalized to 1, as 

— -- = 5(1 - x,)6{l - X2) + a.—-- , (2.5) 

u,X\u,X2 u,X\(JiX2 
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The term of order ag is singular when xi and X2 approach 1. In fact, in this hmit 
the emitted gluon is soft, and one expects a I/X3 singularity. Furthermore, virtual 
corrections, concentrated in the region xi = 1 and X2 = 1 are also infrared- divergent. 
Since our differential cross section is normalized to 1, it is easy to write down a for- 
mula for the average of some physical quantity G{xi,X2) without having to explicitly 
introduce the virtual corrections. We simply write 

{G{xu X2)) = G{1, l) + as dxi dx2 ^j— 7— (^(0:1, X2) - G{1, 1)) . (2.6) 

J (ix\dx2 

The integral is extended in the appropriate Dalitz region. The subtraction term 
under the integral sign embodies the effects of the virtual corrections, and of the 
normalization to the total cross section. Thus, if we choose G{xi,X2) = 1 we clearly 
get 1, which is our normalization. Appropriate formulae for da^^^ / dxidx2 can be 
found for example in ref . . Formula (|2.6|) can be easily implemented in a computer 



program. We computed the following quantities 



where 



f _ _ da^^^ 

Cl = dxi dx2 {xi - 1) , (2.7) 

J axidx2 

f _ _ da^^^ 

di = dxidx2 ,_ ,_ (xia;2-l) , (2.8) 
J dxidx2 

ei = [ dxidx2-ri — 7^ (xi X2 cut(xi, X2) — 1) , (2.9) 
J dxidx2 

cut(xi, X2) = 9 (max(xi, X2) — X3) . (2-10) 



The theta function in the last equation corresponds to the requirement that in doubly 
tagged events the b and the b are in opposite hemispheres with respect to the thrust 
axis. We have 

(x) = 1 + asCi, {X1X2) = I + asdi, (xiX2 cut(xi, 0:2)) = 1 + «sei. (2-11) 
We find, for rrih = 5 GeV, 

Cl = -1.142, di = -2.185, ei = -2.223. (2.12) 



Observe that the coefficients we found are quite large, since one would have naively 
expected them to be of the order of l/vr. This is because large logarithms log E/ nib 
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arise in the coefficients, thus making the perturbative expansion for the above quan- 
tities quite unrehable. The quantities r and r', however, have small perturbative 
coefficients 

,2 

r = ' r, ' ' = as{di — 2ci) = 0.099 , 



(a;ia;2cut(a;i,a;2)) - (a;)' 



a^(ei - 2ci 



0.061a,. 



(2.13) 



Observe that the coefficients are even below the expected magnitude of l/vr. A 
better illustration of the presence of large logarithms, and of their cancellation in 
the quantities r and r' is given in table |I|, where the coefficients Ci, di and ei have 
been computed for m = 10, 1, and 0.1 GeV. While ci, di and ei grow as the mass 
decreases, r and r' approach a finite limit. This is due to the fact that the large 



m 


Cl 


di 


ei 


r 


r' 


10 GeV 


-0.773 


-1.460 


-1.486 


0.086 a. 


0.060 a, 


1 GeV 


-2.053 


-4.000 


-4.045 


0.106 a, 


0.061a, 


0.1 GeV 


-3.348 


-6.590 


-6.636 


0.106 a, 


0.060 a, 



Table 1: Mass depencence of the coefficients Ci, di and ei. 



logarithms present in ci, di and Ci cancel in the linear combinations appearing in 
r. It is easy to prove that this cancellation must occur to all orders in perturbation 
theory. In fact, according to the factorization theorem, the double inclusive cross 
section bb production, in the limit of E ^ rrih, can be written as 

, = / , '^^ D{zi) D{z2) 6{yi zi - xi) 5{y2 Z2 - X2) dyi dzi dy2 dz2 , (2.14) 

dxi dx2 J dyi dy2 

where (in the limit E/mi,^oo) a has a perturbative expansion in a, with finite 
coefficients (observe that the distinction between the momentum- or energy-defined 
Feynman x becomes irrelevant in the limit we are considering here). The divergent 
terms are all absorbed in the fragmentation functions D[z). We then have 

(x) = / dxidx2 xi = ( [ dzzD{z)) [ dyidy2 yi , , (2.15) 

J dxidx2 \J ) J dyidy2 

{X1X2) = J dxidx2 X1X2 ^^^^ = i^j dzzD{z)J J dyidy2 ym ^y^^y ■ (2-16) 
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In the ratio (xi X2) / (x)^ (and therefore in r) the integral containing D cancels. Thus, 
the perturbative coefficients of r are finite in the limit i^/m^— s>oo. Observe that in 
the derivation of eq. (|2.15| ) we have assumed the relation 



D{z)dz = l, (2.17) 

which is appropriate when we can neglect the secondary production of bb pairs via 
gluon splitting. 

Since the cancellation of large logarithms takes place order by order in perturba- 
tion theory, it would be wrong to include incomplete higher-order corrections in the 
expression for r. Thus, for example, if we write 

r = ; = (di - 2ci)as + cWSci - 2di)a. H , (2.18) 

1 -|- CiOg 

the coefficient in front of the term of order is very large when compared to the 
coefficient in front of the term of order as- This is due to the fact that there are 
missing terms in the (9(a^). If a complete 0{a'j.) calculation is performed, these 
large terms cancel. Although a calculation of the next-to-leading term is in principle 
possible, it is not yet available. A calculation of the process of bbgg production has 
been given in ref. [^, but the virtual corrections for the bbg process in the case of 
massive quarks are not available. In order to give an estimate of the next-to-leading 
effects, and to verify the cancellation of the large terms in the 0{a^) coefficients, we 
have computed the 0{a1) contribution to r coming from the real emission process 
bbgg alone. We would like to emphasise that these results can be considered at most 
as an estimate of the order of magnitude of the C(a^) coefficient, since the virtual 
terms are missing. Appropriate subtraction terms have been included in order to 
regulate the soft singularities. We get 

C2 = -2.845 , d2 = -4.595 , 63 = -4.499 . (2.19) 

We have 

r = {di- 2ci)a, + (da - 2cidi - 2c2 + 3cl)a1 + O(a^) 

= 0.099a, + 0.017a^ + C(af) , 
r' = (ei - 2ci)a, + (es - 2ciei - 2c2 + 3c?)a^ + C>(a^) 

= 0.061a, + 0.026a^ + O(a^), (2.20) 
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which shows a well-behaved perturbative expansion. Notice that if we do not include 
the C2, (^2 and 62 terms, the coefficients of the term of order are an order of 
magnitude larger than the coefficient of the term of order a^, which again shows the 
importance of including all enhanced terms. 

The correlation we computed is a small effect, of the order of 1%. In fact, we 
have to decide what is the scale at which as ought to be evaluated in eqs. (|2.13|) . In 
view of our discussion on the cancellation of large logarithms, the correlation should 
be dominated by large momenta, and therefore the appropriate scale should be of 
the order of E. In fact, experience with jet physics suggests the use of a somewhat 
smaller value. For example, the range 0.12 < < 0.16 gives a value of r' between 
0.7 and 1 %. 



3 Double inclusive cross section 



Up to now, we have assumed in our discussion that the efficiency is linear in the B 
momentum. Even in the more realistic case, in which the efficiency is a more compli- 
cated function of the kinematic variables, it is possible to compute the inclusive cross 
section for the production of a hh pair, provided one also knows the B fragmentation 
function, which is to some extent measured at LEP. The appropriate formula is given 
in eq. ( |2.14| ). The expression for the short- distance cross section a up to the order 
is 



da 



5(l-yi)5(l-l/2) + 



2a. 



dyidy2 ' Stt (1 - ?/i)(1 - 2/2) 

The + distribution sign specifies the way that the singularities at yi 
should be treated. For any smooth function of yi, y2 we define 



yi+y2>i 



dyi dy2 



yi)(l -1/2) 



dy\ dy2 



y\ + vl 



Iyi+y2>i (1 - ?/i)(l - ^2) 

Eqs. ( p.l4| ) and (|3.1|) consistently gives 

da f da 

d^, = J ^^^d^ 



G{yi,y2) = 
{Giyi,y2)-Gil,y2)-G{y,,l] 



D{xi) . 



(3.1) 
1 and ?/2 = 1 

(3.2) 
-G(l,l)) . 

(3.3) 



The above equation fixes the factorization scheme to be the annihilation scheme as 
defined in ref. [|T^. Observe that, in this case, the choice of scheme, i.e. eq. 
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fixed unambiguously the result without the need of computing explicitly the virtual 
corrections. In fact, the most general formula for the short distance cross section is 
obtain by adding to eq. (|3.1| ) terms of the form 



^ [5(1 - y,)f{y2) + /(yi)5(l - y2) + g5{l - yi)5{l - y^)] , (3.4) 

TT 

where / is a generic function (in general a distribution) of one variable and g is a 
number. However, in order for eq. (|3.3| ) to be respected, these terms must be absent. 
From eqs. ( ^.14| ), ( p.l| ) and (^) we immediately derive the 0{as) formula 

= D{x^)D{x2) + ^ /'rfi/A%i + y2-l) ^'^.f , X 
dxidx2 37r 7o (1 - ?/i)(l - ^2) 

D (^] D f ^"j ^ e{yi - xi) e{y2 - X2) - D{x,)D (^] - e{y2 - X2) 

\yij \y2j ym \y2j 1/2 

-D(^]D{x2)-e{y^-x^) + D{x,)D{x2) . (3.5) 

\yij yi 

As an illustration, we plot in fig. ^ the double inclusive cross section da/dxidx2 as 
a function of Xi for several values of X2- We use the Peterson parametrization of the 
fragmentation function 

z(l - 

m = N. ,'^ (3.6) 
((1 — 2;)^ + ze)^ 

where is fixed by the condition / D{x)dx = 1. We used the value e = 0.04, which 
gives (x) = 0.70, and = 0.12. The positive momentum correlation is quite visible 
in the figure. As X2 increases, the peak of the distribution in xi also moves towards 
larger values. Using the above formula, we can again compute r. The result should 
not depend upon the choice of the fragmentation function. In this case, the quantity 
(x) does not receive corrections of order a^. We get 

(x)= 0.7036, (X1X2) = 0.4950 + 0.0525 a,, r = 0.1061a,, (3.7) 

which is close to the value previously obtained r = 0.099. In fact, the difference is 
due to mass effects, as can be seen from table |I]. 

There are several ways in which the above leading order calculation could be 
improved in principle. First of all, in the region of xi, X2 near 1, large logarithms of 
1—x arise, and they could be resummed to all order in perturbation theory. It would 
also be desirable to include next-to-leading corrections to the partonic cross section. 
At this moment, these corrections have not yet been computed. 
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Figure 2: Double inclusive cross section da/dxi dx2, plotted as a function 
of xi for several values of X2- 

4 Conclusion 

Monte Carlo programs implement to some extent some of the QCD dynamics 
discussed in the present work. We therefore expect that, in general, at the parton 
level, they should give results consistent with the QCD calculation reported here, 
at least for quantities like r, in which large logarithmic effects cancel. Observe that 
the computation of r must be performed either consistently at some fixed order, or 
including all logarithmically enhanced terms. Failure to do so may expose large, 
uncancelled logarithmic terms, and it may therefore lead to a wrong result. Shower 
Monte Carlo programs include all logarithmically enhanced terms. On the other 
hand, matrix element Monte Carlo programs use a truncated pcrturbative expansion. 
They can be used to compute r only if each fixed-order contribution can be isolated, 
and all terms of order higher than the accuracy of the Monte Carlo can be thrown 
away. 

The QCD factorization theorem guarantees that in quantities like r large logarith- 
mic effects, as well as hadronization effects, should cancel. This statement is valid 
in a leading-twist sense. Thus, it is possible that power-suppressed corrections to r 
are present. A correction of the order of l/Q would be a 1% effect, comparable to 
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the 0{as) corrections. It is difficult to say anything about the importance of power 
corrections, since our knowledge of the hadronization mechanism is quite poor. One 
has therefore to rely upon hadronization models, which are commonly implemented 
in Monte Carlo programs. The formalism presented in this paper can also be applied 
to the computation of other processes, such as, for example, the correlation in the 
inclusive production of strange mesons in the large-x region. These processes may 
provide a testing ground for the theoretical predictions, and may help us gain confi- 
dence on the reliability of the hadronization models implemented in the Monte Carlo 
programs. 
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